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An exact method for direct calculation of the Jost function and Jost solutions for a repulsive 
singular potential is presented. Within this method the Schrodinger equation is replaced by an 
equivalent system of linear first-order differential equations, which after complex rotation, can eas- 
ily be solved numerically. The Jost function can be obtained to any desired accuracy for all complex 
momenta of physical interest, including the spectral points corresponding to bound and resonant 
states. The method can also be used in the complex angular-momentum plane to calculate the 
Regge trajectories. The effectiveness of the method is demonstrated using the Lennard-Jones (12,6) 
potential. The spectral properties of the realistic inter-atomic 4 He- 4 He potentials HFDHE2 and 
HFD-B of Aziz and collaborators are also investigated. 
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I. INTRODUCTION 

A new method for locating potential resonances and Regge trajectories, based on direct calculation of the Jost 
function in the complex fc-plane, has recently been developed |l|-|4|. Within this method, the bound, resonant, and 
scattering states can be found by calculating the Jost solutions and the Jost function on the appropriate domain of 
^C} ' the fc-plane. The bound and resonant state energies, for example, can be found by locating the zeros of the Jost 
£N) , function on the positive imaginary axis and in the fourth quadrant respectively. At the same time, as a by-product 
of the Jost function calculation, one gets the physical wave function that has the correct asymptotic behavior. 

t-h ; 

The method, in the form developed in Refs. Jl]-Q, cannot be directly applied to potentials which are more singular 
than 1/r 2 at the origin. A significant number of practical problems, however, where the Jost function could be very 
useful, involves such potentials. For example, inter-atomic and inter-molecular forces at short distances are strongly 
repulsive due to the overlap of the electron clouds, and thus they are usually represented by repulsive singular poten- 
tials such as the Lennard-Jones(12,6) one which has ~ 1/r 12 behavior. It is, therefore, desirable to extend the Jost 
' function method for potentials of this kind. 
3 ' 

It is well-known that attractive singular potentials do not admit physically meaningful solutions with the usual 
boundary conditions || . All solutions of the Schrodinger equation with such potentials vanish at the origin and there 
is no apparent way to determine the arbitrary phase factor between them. In contrast, repulsive singular potentials 
do not pose any problem regarding mathematical uniqueness or physical interpretation. However, the integration of 
the Schrodinger equation as well as of the relevant equations for the Jost function, have inherent difficulties resulting 
from the fact that the singularity of the potential makes the r — an irregular singular point of the equation. In 
particular, the regular solution cannot be defined by universal boundary conditions independent of the potential. 

This drawback, stemming from the extremely strong repulsion near the origin, can fortunately be tackled by the 
WKB approximation which provides the correct radial behavior of the wave function in the neighborhood of the point 
r = ^jj. Therefore, starting with the WKB boundary conditions, one can find the regular solution by integrating 
the Schrodinger equation from r — to some intermediate point r = r- mt . Then the equations for the Jost function 
can be integrated, from r- lnt outwards as for the nonsingular potential using the boundary conditions at r = r- mt which 
are expressed in terms of the regular solution and its first derivative. 

The paper is organized as follows: In Sec. II our formalism is presented and is tested in Sec. Ill using an 
example known in the literature; in Sec. IV the method is applied to realistic interaction between helium atoms. Our 
conclusions are drawn in Sec. V. 
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II. THEORY 



A. Basic equations and definitions 

There are three different types of physical problems associated with the Schrodinger equation (2/i/h 2 = 1), 

[8 2 r + k 2 ~ £(£ + l)/r 2 ] u t {k, r) = V{r)u t {k, r) , (1) 

namely, bound, scattering, and resonant state problems. They differ in the boundary conditions imposed on the wave 
function at large distances. Alternatively, a solution can be prescribed by the boundary conditions at the origin. In 
the case of regular potentials obeying the condition 

lim r 2 V(r) = , (2) 
r— >0 

the solution </>e(k,r) which vanishes near r = exactly like the Riccati-Bessel function, 

lim </)£(k,r)/ji(kr) = 1 , (3) 

r— >0 

is called the regular solution. Since all physical solutions are regular at the origin, they differ from (f>£(k,r) only by a 
normalization constant. Therefore, if the function (f>i(k,r) can be calculated at all real and complex momenta k, one 
can have, in principle, all solutions of physical interest in a most general form. For example, the calculation of the 
scattering solutions on the real fc-axis is simply a matter of finding the proper normalization for </>^(fc,r) because the 
regular solution has the correct behavior at large r for any k > 0. In contrast, to find the bound and resonant states 
where the k is complex, one must ensure that the function 4>i(k, r) has the proper physical asymptotic behavior which 
exist only at certain points on the fc-plane. These spectral points can be found by many different ways, but perhaps 
the most convenient way to find them is by locating the zeros of the Jost function. 

For any complex k, the regular solution at large distances can be expressed as a linear combination of the Riccati- 
Hankel functions h^\kr), 

Mk,r) — ► \ \h { e +) (kr)f*(k*) + h { -\kr)f e (k)} , (4) 

where the r-independent but momentum-dependent coefficient fe(k) is the Jost function. From the asymptotic form 
(||) it is clear that the zeros of fe(k) on the positive imaginary axis of the complex fc -plane correspond to bound states 
while those in the fourth quadrant to resonances. 

In order to find fi(k), we look for the regular solution on the whole interval [0, oo) in the form 

<f>i(k,r) = ^ [h i +\kr)F^ +) (k,r) + h ( f\kr)F^\k 1 r)\ , (5) 

where the new unknown functions F^\k, r) are subjected to the additional condition 

h ( e +) (kr)d r F^ + \k,r) + h ( f\kr)d r F^\k,r) = . (6) 
Eq. (|l|) is then transformed into an equivalent system of first order equations 

d r F^\k,r)^± H \^ r) V(r) [^ +) (fcr)F £ (+) (fc,r) + (kr)F^\k,r)} . (7) 

In Refs. HH it was shown that at large distances F^ \k, r) coincides with the Jost function, 

lim F { -\k,r) =h{k) , (8) 

r— >oo 

but this limit only exists when 

Im kr > . (9) 

If r is real, the condition (^|) is only satisfied for bound and scattering states but not for resonances. To calculate 
fe(k) we, therefore, make a complex rotation of the coordinate in Eqs. (0), in the first quadrant 

r = xexp(i6), x>0, 0<6»<|, (10) 



with a sufficiently large 6 (see Refs. |l|-|4|] for more details). Such a rotation is only possible if the potential is an 
analytic function of r and tends to zero when x — > oo for the chosen angle 9. 
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B. Boundary conditions 



In the case of regular potentials the boundary conditions for Eqs. (Q) are very simple, 

F] ±) (fc,0) = l. (11) 
They follow immediately from (||), (||), and the fact that 

i [h ( +\kr) + h<f ) (kr)\ =j t {kr) . 

Going over to singular potentials, Eq. pi) does not hold anymore. Due to the extremely strong repulsion, the regular 
solution vanishes much faster than ji{kr) when r — > 0. In fact, it vanishes exponentially |7j] and therefore the conditions 
( pi] ) must be modified accordingly. In order to find the exact behavior of the regular solution near the origin we apply 
the familiar semi-classical WKB method. Though the strong repulsion makes things rather complicated, it has the 
advantage that the criterion of the applicability of the WKB approximation is satisfied when r — > 0. Indeed, the 
WKB method works well when the local wavelength A varies slowly, i.e. 



\d\/dr\ < 1 



It can be shown [|| that this derivative is given by 

|dA/drj = • 



^V-V(r)]- 3/2 
dr 



Assuming that V(r) approaches its singularity near r = monotonically, we can find an r n 
momentum in (|l3|) is negligible, i.e. we may write 



(12) 



(13) 



that for all r < r m i„ the 



|dA/dr 



1 

2 



dV(r) 



dr 



[V(r)} 



-3/2 



(14) 



When r — * 0, the right hand side of Eq. ( |i~4| ) for usual singular potentials tends to zero. For example, if 



V(r) — ► s /r B 

r— »0 



the condition ([L2|) is always satisfied for n > 2, 

|dA/dr| - 



in-l 

nr 2 
2^5 



, 



if n > 2 



Therefore, assuming that the necessary condition ( |i~2| ) is fulfilled and choosing a small enough r m ini w e can express 
the regular solution on the interval [0, r m \ n ] using the WKB approximation (see, for example, Ref. S), viz. 



<t»t(k,r) 



exp 



i / p{p)dp 



r e [0, r n 



where the classical momentum p(r) is defined by 

p(r) = 



/c 2 -^(r)-(£+-) 2 /r 2 



(15) 



(16) 



and the upper limit a in the integral is an arbitrary value a > r m i„ . Usually a is taken to be the inner turning point 
Hi, but it is obvious from Eq. ([l5]) that an additional integration from a to the turning point can only change the 
overall normalization of the solution which is not our concern at the moment. Thus, Eq. (|l^) together with the 
derivative 



d r (pi(k,r) 



dV(r) 
dr 



2[£ + 



1 



4[p(r)] 5 /2 



exp 



p(p)dp 



r e (0, r min ] 



(17) 
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can be used as boundary conditions for the regular solution of the Schrodinger equation at any point in the interval 
(0,r m ; n ]. To obtain the corresponding boundary conditions for the functions F e (k,r), we need to express them in 
terms of <fie(k,r) and d r <fie(k,r). For this we can use Eq. (||) together with relation 



dAi{k,r) = - \Fl + \k,r)d r hY>{kr)+F ( t '{k,r)d r h\ >(kr) 



,(+)/ 



,(-)/ 



which follows from (||). From (||) and (|T^) we find that 



F < (±) (fc,r) = ±i [<j> l {k,r)d r hf\kr) - h!f\kr)d r <i>i{k,r) 



(18) 



(19) 



which is valid for any r G [0, oo). Therefore Eqs. ( [l9| ) taken at some point r < r m ; n with cj>t(k, r) and d r <pi(k, r) given 
by ( |l5| ) and (|l7|), provide us the boundary conditions, required in Eqs. f^), for singular potentials. It can easily be 
checked (by using ji{kr) for the regular solution near r = 0) that Eq. (|H) gives the correct boundary conditions for 
regular potentials as well, 

Alternatively to impose the boundary conditions on the functions F^\k,r) near the origin, one can simply solve 

the Schrodinger equation from a small r up to some intermediate point b where, using pi), the F^\k,b) can be 
obtained and propagated further on by integrating equations (^). 



C. Integration path 

The use of more complicated boundary conditions at r — does not change the condition (Q) for the existence of 
the limit (^|). Indeed, in deriving this condition we used only the behavior of the potential and the Riccati-Hankel 
functions at large distances plH]. Therefore, the Jost function for a singular potential can also be calculated by 
evaluating the function Fj \k,r) at a large r. When we are dealing with resonances, i.e. working in the fourth 
quadrant of the fc-plane, we need to integrate Eqs. (Q) along the turned ray (|To|). 



As can be seen from the WKB boundary conditions (|l£j), the use of a complex r near the origin, makes </>^(A;,r) 
oscillatory from the outset. Although this does not formally cause any problem, in numerical calculations such 
oscillations may reduce the accuracy. To avoid this we solve Eqs. (^) from a small r m j n to some intermediate point b 
along the real axis and then perform the complex rotation, 

r = b + x exp(i0) , xe[0,oo), < 6 < | , (20) 

as is shown in Fig. [l]. Therefore, on the interval [r m j n ,6] we can use Eqs. (0) as they are, while beyond the point 
r = b these equations are transformed to 

d x Ff\k, b + xe w ) = ± eteh P(^+ kxet9 ) v{b + x ^ Yf\kb + kxe ie )F^ + \k, b + xe ie ) 

(21) 

+ hi" \kb + kxe l8 )F^ ^k^b + xe 10 ) . 



Though the complex transformation (|20|) is different from (JlOj), the proof of the existence of the limit (|8j) given in 
the Appendix A. 2 of Ref. Q remains applicable here. Indeed, that proof was based on the fact that for Imfcr > 
the Riccati-Hankel function h^\kr) decays exponentially at large |r|, and thus the derivative d r F^(k,r) vanishes 
there and the function F^~\k, r) becomes a constant. Under the transformation j20| ) the asymptotic behavior of the 
Riccati-Hankel function, 

h ( +\kr) — ► -iexp[i(kr-l-K/2)] , (22) 

T — >00 

has only an additional r-independent phase factor exp(ikb) which does not affect the proof. 

From the above, it is clear that we can identify the Jost function fi(k) as the value of Fg \k,b + xe %9 ) at a 
sufficiently large x beyond which this function is practically constant. In the bound and scattering state domain, 
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where Imfc > 0, one can choose any rotation angle 9 allowed by the potential, including 9 = 0. In the resonance 
domain, however, where 

k = \k\ exp(— i</?) , if> > , 

the rotation angle 9 must be greater or equal to <p. If the condition 9 > ip is fulfilled, the value of the limit (||) does not 
depend on the choice of 9. This provides us with a reliable way to check the stability and accuracy of the calculations 
by comparing the results for fi(k) obtained with two different values of 9. 

From Eq. it is clear that the angular momentum appears only in the phase factor of the asymptotic behavior 
of the Riccati-Hankel functions and hence of the regular solution. Therefore, the use of any complex £ cannot change 
the domain of the /c-plane where the limit (||) exists. This means that the Jost function can be calculated, for any 
complex angular momentum, using the same equations. Moreover, when looking for the Regge poles in the ^-plane, 
the complex rotation is not necessary because these poles correspond to real energies. Locating Regge poles as zeros 
of the Jost function in the complex £-plane is easier than by calculating them via the S'-matrix using three integration 
paths (in the r-plane) as suggested in Ref. ||. 



D. Jost solutions 



By storing the values of iy (k, r) on the integration grid one can also obtain the regular solution in the form (|^) 
on the interval [r m ; n , r max ]. It is noted that the use of the Riccati-Hankel functions in (|^) guarantees the correct (in 
fact exact) asymptotic behavior of the wave function. 



The regular solution thus obtained consists of two terms: 



i/4 +) (fcr)^ (+) (fc,r) 



■exp[+i(*cr-*7r/2)] //(*:*) > 



h { -\kr)F^'(k,r) 



(-), 



+ -exp[-i(kr-£Tr/2)]ft(k) 



Asymptotically they behave like ~ e ±lfcr and thus at long distances they are proportional to the commonly used Jost 
solutions f^\k, r) for which 



ff } (fc,r) 



hf\kr) 



(23) 



In practice, the Jost solutions can be calculated, via (g), by integrating Eqs. (Q) inwards from a sufficiently large r n 
with the boundary conditions 



F^ +) (k r ) 

± £ V ' 'max ) 

Fg (&, r max ) 





' 2 " 








for fj +) (fc,r) , 



(^3 ''"max) 





' " 




2 



for 



A \k,r) 



which obviously comply with the definition (E3J). The advantage of such an approach is that at large r all the 



oscillations of f^ (fc, r) are described exactly by the Riccati-Hankel functions while the functions Fg ± \k,'i 



arc 



smooth. 



III. LENNARD-JONES POTENTIAL 



In order to evaluate the accuracy and efficiency of our method we apply it to the Lennard- Jones potential 

12 / ,\ 6" 



V(r) = D 



r 



(24) 
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which is well-known in atomic and molecular physics. Combined with a rotational barrier, this potential supports 
narrow as well as broad resonant states (see, for example, Ref. §). To locate them, any method employed must be 
pushed to the extreme, thus exhibiting its advantages and drawbacks. 

To be able to compare our results with other calculations, we chose the parameters in ( p4| ) to be the same as those 
used in Refs. |§Jl0||, namely, d = 3.56 Aand with D varying from 5cm -1 to 60 cm -1 . The choice D — 60 cm -1 together 
with the conversion factor h 2 /2fi = 8.7802375 cm 1 A 2 (which was used for all values of D) approximately represents 
the interaction between the Ar atom and the H2 molecule [||. 

In Tables | and |l| the energies and widths of the first resonant states in the partial wave i = 8 are presented for 
different values of D. The results obtained with three other methods described in Rcfs. [^],[l0| are also given. The digits 
shown there are stable under changes of the rotation angle and thus they indicate the accuracy achieved. The third 
column of these tables, contains the results obtained in Ref. || using a Complex Rotation (CR) method which in some 
aspects is similar to ours. The authors of that reference perform the rotation directly in the Schrodinger equation and 
integrate it from r = outwards and from a large r max inwards. At the origin they use the WKB boundary conditions 
and at r max they start from the Siegert spherical wave. In other words, the wave function is calculated using physical 
boundary conditions. In such an approach a resonance corresponds to a complex energy which matches the inward 
and outward integration. As indicated in Ref. this method fails for broad resonances due to instability in the 
outward integration. In the fourth column the results obtained in Ref. Q using the Quantum Time Delay (QTD) 
method are cited. This method is expected to be reliable for narrow resonances but its applicability to broad states 
is questionable. Finally, in the last column of Table | and |l| we give the results obtained in Ref. using the Finite 
Range Scattering Wave (FRSW) method. The main idea of this method is based on the fact that while the scattering 
wave function cannot be expanded properly by a finite number of square integrable functions on an infinite range, it 
is possible to do so for a finite range. 

The test calculations show that our method works well, especially for narrow resonances. Broad resonances can 
also be located. In contrast to the CR-method of Ref. 0, which was unstable for broad resonances corresponding 
to D < 35 cm -1 , we succeeded even in the case of D = 5 cm -1 which generates an extremely broad state (its width 
is greater than the resonance energy by a factor of 2). Our results for small values of D, reproduce well the curve 
depicted in Fig. 3 of Ref. ||] which was produced semi-classically. The greater stability of the Jost function method 
as compared to the CR-method of Ref. || can be attributed to the use of the ansatz (||) for the regular solution. 
The Riccati-Hankel functions, explicitly extracted there, describe correctly all oscillations at large distances with the 
remaining functions Fg being smooth. Another reason for this stability is the use of the deformed integration path 
shown on Fig. El which enables us to avoid fast oscillations at short distances. 



IV. AZIZ POTENTIALS 



The model potential considered in the previous section, though of typical form for inter-molecular interactions, does 
not describe any real physical system. To give a more practical example, we apply our formalism to study the inter- 
action between two 4 He atoms. This interaction is of interest in the Bosc-Einstein condensation and super-fluidity 
of helium at extremely low temperatures. It is known that two helium atoms form a dimmer molecule with binding 
energy of ~ lmK, but, to the best of our knowledge, the possibility of forming dimmer resonances has not been 
investigated yet. 



The search for a realistic 4 He- 4 He potential is a long-standing problem in molecular physics. The earliest successful 
potential of the Lennard-Jones(12-6) form was fitted just to reproduce the second virial coefficient. Later on some 
other characteristics of helium gas, such as viscosity, were included into the fitting (for a more detail review see Refs. 
flTiyni ). Nowadays, the potentials suggested by Aziz and co-workers are considered as be realistic. Therefore, in this 
section, we apply our method using two versions of these potentials, namely, the HFDHE2 and the HFD-B |l2] ] 
potentials. They can both be described using the same analytical form 




^exp(-aC-/?C 2 ) 



-(B/c-iy 



if C < B 
if C > B 



C10 



no 



(25) 
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but with different choices of the parameters (see Table [II). The only principal difference in the functional form 
between them is the absence of the Gaussian term ((3 = 0) in the HFDHE2 potential. 



Formally, the HFDHE2 and HFD-B are regular potentials since the presence of the cut-off function F(() in 
makes them finite at r = 0, 

V(r) —> eA . (26) 

However, the product eA is very large (~ 10 6 ) as compared with the values of the potential in the attractive region. 
This causes numerical instabilities when one tries to solve the Schrodinger equation using methods designed for regular 
potentials. To avoid this difficulty, we notice that like in the case of singular potentials the fast growth of the repulsion 
near the origin allows the use of the WKB boundary conditions near r = 0. Indeed, the derivative of the potential in 
the vicinity of this point, 

^ - , (27) 

dr r->o r m 

is of the same order of magnitude as V , which makes the derivative of the local wavelength ( |l4| ) very small because 
of the large A, 

\d\/dr\ -j °== , (28) 

2r m J-^eA 



where the conversion factor h 2 /2fi — 12.12 KA 2 , corresponding to the choice of the units in Table III. should be used. 



With the parameters given, formula Q28|) gives 0.003 and 0.004 for the potentials HFDHE2 and HFD-B respectively. 



These values of |dA/dr| are small enough to comply with (12) and allow the use of WKB boundary conditions. We 



can, therefore, apply the method described in the preceding sections, to the potentials HFDHE2 and HFD-B as if 
they were singular potentials. 

To begin with, we tested the ability of our method to deal with thi s kind of potentials by calculating the dimmer 



binding energy. The results of these calculations are given in Table IV where, for comparison, we also cite the binding 
energies obtained in several earlier works. It is seen that the potentials HFDHE2 and HFD-B support a dimmer 
bound state at energies which differ by a factor of 2. 

A question then arises whether these potentials generate also quite different distribution of resonances which would 
result in different on and off the energy shell characteristics of the scattering amplitude. To study this we located 
several zeros of the Jost function in the momentum as well as in the ^-plane (Regge poles) for both potentials. Due 
to the absence of a potential barrier there are no resonances in the S-wave (at least with a reasonably small width) . 
They appear, however, at higher partial waves, starting from £ — 1. The energies and widths of several such reso- 
nances are given in Table They are the lowest resonant states in each partial wave as they belong to the same 
Regge trajectory which starts from the ground state. The trajectories for the potentials HFDHE2 and HFD-B are 
practically indistinguishable and are shown in Fig. @ by a single curve. Few points of this curve which correspond to 



resonances, are also given in Table |VI[ It is seen that, to all practical purposes, the position of the Regge poles are 
the same. 

As can be seen in Table |v|, in each partial wave the potential generates a broad resonance which covers the whole 
low-energy region. This, together with the fact that the bound state pole of the amplitude is very close to k = 0, 
implies that the cross-section at energies ~ 10 °K (~ 10~ 3 eV) must be quite large. 



V. CONCLUSIONS 



We presented an exact method for calculating the Jost solutions and the Jost function for singular potentials, for 
real or complex momenta of physical interest. We demonstrated in the examples considered, the suggested method is 
sufficiently stable and effective not only in the case of true singular potentials but also when a potential has strong, 
though finite, repulsion at short distances. 
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The method is based on simple differential equations of the first order, which can be easily solved numerically. 
Thus, the spectrum generated by any given potential can be thoroughly investigated. At the same time, physi- 
cal wave function can be obtained having the correct asymptotic behavior. When the potential has a Coulomb 
tail one can simply replace the Riccati-Hankel functions in the relevant equations by their Coulomb analogous, 
Hj^fy, kr) = Fi(r), kr) =p iGg(r], kr) Q. In the case of a non-central potential the Jost function as well as the differ- 
ential equations assume a matrix form with somewhat more complicated, but still tractable boundary conditions at 
r = 0|. 

The method is also applicable when the angular momentum is complex. This enables us to locate Regge trajectories 
as well. This could be useful, for example, in molecular scattering problems where the partial wave series in many 
cases converges slowly [fL7| . This slow convergence can be overcome by allowing the angular momentum to become 
complex valued which allows the use of the Watson transformation. However, such a procedure requires the knowledge 
of the positions of the Regge poles. 
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Ref. 


This work 


CR § 


QTD | 


FRSW M 

E rea (cm" ) 


D (cm" 1 ) 


E TCS (cm -1 ) 


S res (cm 1 ) 


_B ros (cm 1 ) 


5 


30 








10 


27 








15 


25.5 








20 


24.5 








25 


22.90 








30 


21.193 








oo 


i q a 


1 Q /1/1Q 


1 Q ^70 




40 


17.6478 


17.647 


17.619 


17.617 


45 


15.7768 


15.777 


15.769 


15.769 


50 


13.81980 


13.820 


13.819 


13.818 


55 


11.744242 


11.744 


11.744 


11.743 


60 


9.4943275 


9.494 


9.494 


9.4934 



TABLE I. Energies of the lowest resonances, in the £ — 8 partial wave, for the Lennard-Jones potential with different D. 



Ref. 


This work 


CR § 


QTD | 


FRSW M 


D (cm" 1 ) 


r res (cm -1 ) 


r rcs (cm 1 ) 


r rcs (cm 1 ) 


r rcs (cnr^) 


5 


60 








10 


42 








15 


29.4 








20 


24.8 








25 


18.63 








30 


13.70 








35 


9.724 


9.727 


10.228 




40 


6.533 


6.536 


6.604 


6.603 


45 


4.039 


4.039 


3.992 


3.990 


50 


2.1833 


2.183 


2.143 


2.142 


55 


0.93915 


0.939 


0.926 


0.926 


60 


0.264474 


0.264 


0.263 


0.263 



TABLE II. Widths of the lowest resonances, in the £ = 8 partial wave, for the Lennard-Jones potential with different D. 



parameter 


HFDHE2 


HFD-B 


e(K) 


10.8 


10.948 


r m (A) 


2.9673 


2.963 


A 


544850.4 


184431.01 


a 


13.353384 


10.43329537 


P 





2.27965105 


Ce 


1.3732412 


1.36745214 


Ce 


0.4253785 


0.42123807 


Cio 


0.178100 


0.17473318 


B 


1.241314 


1.4826 



TABLE III. Parameters of the two versions of the Aziz 4 He- 4 He potential. 
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4 He- 4 He binding energy (mK) 



Ref. 



HFDHE2 



HFD-B 



This work 
P 
[12 
[14 
[15 
[16 



0.8301249029 
0.8301 

0.830 
0.829 
0.8299 



1.6854110471 

1.6854 

1.684 



TABLE IV. Binding energies of He2 di-atomic molecule for the two versions of the Aziz potential. 



HFDHE2 



E(K) 



r (K) 



E(K) 



HFD-B 



r(K) 



0.334 
2.164 
5.930 
11.954 
20.478 



1.822 
6.825 
15.195 
27.196 
42.949 



0.339 
2.179 
5.963 
12.005 
20.518 



1.795 
6.774 
15.098 
27.000 
42.603 



TABLE V. Energies and widths of the lowest resonant states generated by the two versions of the Aziz potential in several 
partial waves. 



HFDHE2 


HFD-B 


E (K) 


£ 


E(K) 


£ 


0.334 




0.525 -M0.429 


0.339 




0.535 + i0.424 


2.164 




1.408 + i0.778 


2.179 




1.416 + i0.773 


5.930 




2.294 + il.065 


5.963 




2.303 + il.060 


11.954 




3.182 + il.332 


12.005 




3.191 326 


20.478 




4.070 + il.590 


20.518 




4.078 + il.584 



TABLE VI. Regge poles corresponding to resonances generated by the two versions of the Aziz potential. 
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FIG. 1. Deformed contour for integration of the differential equations 
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FIG. 2. Regge trajectory for the HFD-B potential. Filled circles indicate bound and resonant states. 
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